Evolution of Spin-Orbital-Lattice Coupling in the RVO^ Perovskites 
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We introduce a microscopic model which unravels the physical mechanisms responsible for the 
observed phase diagram of the RVO3 perovskites. It reveals a nontrivial interplay between superex- 
change, the orbital-lattice coupling due to the GdFeOa-like rotations of the VOe octahedra, and 
orthorhombic lattice distortions. We find that the lattice strain affects the onset of the magnetic 
and orbital order by partial suppression of orbital fluctuations. The present approach provides also 
a natural explanation of the observed reduction of magnon energies from LaVOs to YVO3. 
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Over the last decade, extensive work on transition 
metal oxides has demonstrated a strong interrelationship 
between spin order (SO) and orbital order (00), often 
compounded by the occurrence of various types of lat- 
tice distortions, resulting in phase behavior of dazzling 
complexity. Recently, however, impressive experimental 
work has produced exceptionally detailed information on 
the phase diagrams of the RMriO^ manganites ll| and 
the i^VOa vanadates (where i?=Lu,Yb,- • •,La) [2|, thus 
providing a unique challenge to the theory and the op- 
portunity to resolve the interplay between the underlying 
microscopic mechanisms. 

The manganite i?Mn03 perovskites exhibit the more 
common behavior, i.e., upon lowering the temperature, 
the 00 appears first, accompanied by a lattice distortion, 
at Too I s-iid thus modifies the conditions for the onset of 
SO at a significantly lower temperature Tn- When the 
ionic radius rn of the B?'^ ions decreases, the Neel tem- 
perature Tjv gets drastically reduced and the 00 tran- 
sition temperature Too is enhanced [ij. By contrast, 
in the BNO3 vanadates the two transitions are close to 
each other, i.e., Tjvi ^ Too, the type of order being G- 
type 00 (G-OO) and C-type antiferromagnetic (C-AF), 
setting in below Too and T^i 0]i respectively [J]. More- 
over, whereas T^vi again gets reduced for decreasing r^, 
Too exhibits a nonmonotonic dependence on r^j 0]. 

These experimental results demonstrate that spins and 
orbitals are intimately coupled in the RYOz vanadates, 
consistent with the recent observation that these com- 
pounds form a unique class characterized by strong or- 
bital fluctuations 0, H 01 which follow from superex- 
change interactions between almost degenerate t2g or- 
bitals and spin-orbit term [HI • Their coupling 
is also visible in spectacular changes of the SO and 00 
under pressure The pressure dependence and ther- 
mal conductivity data [l^ suggest in turn strong orbital- 
lattice coupling 13|. As in t2g systems Jahn- Teller (JT) 



for being involved in the coupling between orbitals and 
the lattice. 

In this Letter we present a first microscopic approach 
to the phase diagram of the KVO^ vanadates using an ex- 
tended spin-orbital model which treats the coupled spin 
and orbital degrees of freedom in the presence of lat- 
tice distortions. We demonstrate that the generic trends 
observed in the phase diagram of i?V03 can be indeed 
explained by the theory, see Fig. [U provided one in- 
cludes explicitly the coupling between the orbitals and 
the lattice distortions that increase with decreasing r^. 

A priori, the decreasing ionic radius rn in the i^VOs 
perovskites triggers the GFOD obtained by alternating 
rotations of the VOg octahedra by an angle d around the 




interactions are expected to be rather weak, the GdFeOa- 
like distortions (GFOD) [li, 15| are the prime candidate 



FIG. 1: (color) The orbital transition Too and Neel Tni tem- 
perature (solid lines) for varying ionic size in BNO3, as ob- 
tained from the present theory for the parameter values given 
in the text, and from experiment (full and empty circles) Q]. 
Dashed lines indicate Too and Tjvi obtained under neglect 
of orbital- lattice coupling [g^s = 0). The inset shows the 
GdFeOa-type distortion, with the rotation angles ■& and </? 
corresponding to the data of YVO3 [l6| . 
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b axis, and by an angle ip around the c axis (see inset 
in Fig. [T|). This results in a decrease of V-O-V bond 
angles (e.g. = tt — 2i? along the c direction), and leads 
to an orthorhombic distortion u = (b — a) /a, where a 
and b are the lattice parameters of the Pbnm structure. 
Although these changes are systematic in the oxides [17j , 
their relation to the evolution of microscopic parameters 
and physical properties such as the onset of 00 and SO 
along the RVO3 series is not yet fully understood. 

The spin-orbital model for RYO3 includes: (i) the su- 
perexchange interaction 0], (ii) the crystal field (CF) 
splitting which follows from the GFOD, (iii) orbital- 
orbital interactions induced by orbital-lattice coupling, 
(iv) orbital-strain coupling. It takes the generic form 
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{ij) 

+ i?,(z?) J^e'^'-S^f + Kfc(^?) E ^tr^+Uu, (1) 

i {ij)± c 

where the parameters {E^, Vat, Vc} depend on the tilting 
angle d. The first term describes the superexchange of 
strength J — At^ /U [here t is the effective {ddTi) hopping 
between t2g orbitals and U the intraorbital Coulomb in- 
teraction] between V^^ ions in the cP configuration with 
S = 1 spins. The dependence of J on the rare earth 
ion R is weak and is neglected in the present theory; 
we adopted J = 202 K for the theoretical curves in Fig. 
[TJ The orbital operators Jij and /Cy follow from virtual 
dIcPj d^d}j charge excitations and depend on Hund's 
exchange parameter Jh/U. Their form depends on the 
(ij)-bond orientation; they are given in Ref. 9] for the 
actual {xyY(yz/zx)^ configuration in cubic vanadates. 
The orbital (pseudospin) operators = \{nyz — nzx)i 
refer to the two active orbitals: yz and zx [1, Q. The 
CF splitting of these two orbitals oc E^ supports C-type 
00 t^, with a modulation vector Q = (7r,7r, 0) in cubic 
notation. The Vat > and Vc > orbital interactions 
are due to the JT and GFOD distortions of the VOg oc- 
tahedra, and like E^ favor C-type 00. Unlike for Vat, 
it may be expected that the dependence of Vc on the 
angle -d is weak, and indeed a constant Vd'd) = 0.26 J re- 
produces a simultaneous onset of SO and 00 in LaVOs 
within the present model 18], i.e.. Too = Tni, see Fig. 
[TJ Finally, Hu describes the coupling of the orbitals to 
the orthorhombic distortion u and is explained below. 

To derive the functional dependence of the microscopic 
parameters {Ez,Vab} on the tilting angle we consid- 
ered the point charge model, and used the structural data 
for i?V03 19]. Due to the GFOD shown in Fig. [H the 
two active yz / zx orbitals experience the CF splitting E^ ■ 
By considering the ionic charges acting on the t2g orbitals 
and analyzing the level splittings, we obtained: 



Ez{^9) ^ Jvz sin^iJcosi?, 



(2) 



the {yz, zx} orbitals. These qualitative trends are valid 
in a range of and for further analysis we adopted a 
representative value — t?/2, similar to the trend in 



titanates 



151 ]. It is expected that the angular dependence 



of the orbital interaction Vat follows the CF term ^ : 

VabW = JVab sin^l^COSl?. 



(3) 



An important term in ([T]), coupfing the orbitals to the 
lattice, is the one involving the orthorhombic strain u, 



Till 



--NK{u~uom\ 



(4) 



where 5 > is the coupling constant, K is the force 
constant, and N is the number of V^~^ ions. In contrast 
to the longitudinal field E^, gu acts as a transverse field, 
i.e., it favors that one of the two linear combinations 
-^{\yz)zL\zx)) is occupied. Since u is a classical variable, 
we may minimize Eq. Q and write the global distortion 
as ul-d; T) = uq{{)) + [g/ K){t^)t, consisting of (i) a pure 
lattice contribution uq{§), and (ii) a contribution due 
the orbital polarization cx (r^) which we determined self- 
consistently. Both uq and (r^) are expected to increase 
with increasing tilting As will be shown below, (r^) 
has only a weak temperature dependence, so we may use 



gca{'&) = gu{d) — Jvg sin^i^cosi^. 



(5) 



Indeed, we shall see below that this strong dependence of 
gcs on -d is not only necessary to reduce Too for large tilt- 
ing angles i?, but is also consistent with the experimental 



data for m(i?) [H, S [IH, [H. Altogether, {v 
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and verified that the xy orbitals are indeed well below 



are the parameters of the spin-orbital model ((T]) 

In the i?Mn03 manganites, mean-field (MF) theory 
with classical, on-site, SO and 00 parameters can be 
used to investigate the phase diagram [23|]. However, this 
approach with on-site order parameters does not suffice 
in the vanadates 2j] when orbital fluctuations stabiliz- 
ing the C-AF phase are present Q — then it becomes 
essential to determine self-consistently the orbital singlet 
correlations (t^ ■ fj) as well. Hence we used a cluster MF 
theory for a bond {ij) along the c axis [2^, with spin 
and orbital MF terms (5"=) and {t'-)g = ^Kt/ - tJ)] 
which follow from interactions with the neighbors in 
all three cubic directions. Apart from the singlet orbital 
correlations {fi ■ fj), the transverse field (r^) was crucial 
to reproduce the phase diagram of Fig. [T] 

The nontrivial character of the phase diagram and the 
underlying spin-orbital coupling in the i?V03 vanadates 
can be fully appreciated by analyzing the variation of 
the microscopic interactions with decreasing angle Q (in- 
creasing tilting 1?) . While the CF splitting and the orbital 
interactions Vab show only moderate increase for decreas- 
ing O, the orbital polarization geff increases rapidly and 
becomes quite large when G < 150° (Fig. Note that 
the increasing JT term Vab supports the superexchange 
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FIG. 2: (color online) Parameters {Ez,Vab, dcs} in Eq. ((T| 
(in units of J) for varying V-O-V bond angle O (in degrees). 
Parameters: — 17, Vab = 22, v„ = 740. 



and stabilizes the G-00, both the increasing CF splitting 
Ez, and the orbital-lattice coupling g^ff compete with it. 

With the present parameters 00 and SO occur simul- 
taneously in LaVOa, and Too = Tjvi - 0.73J (Fig. ^. 
The orbital polarization (r^) ~ 0.03 is here rather weak 
at Tjvi, and is further reduced in the ordered phase when 
the 00 parameter {t^)g grows with decreasing T < Tqo 
(due to finite E^, the orbitals xz/zy are nonequivalent 
and (r^)o = Kr/)| > even for T > Too [11]). In con- 
trast, in SmVOs the 00 occurs first at Too — 0.86 J, and 
the SO follows only at T/vi ~ 0.65J. For these param- 
eters the transverse orbital polarization is considerably 
increased, with (t^) ~ 0.20 at Tni (see Fig. [3]). Note 
that the polarization (r^) does not change at T ~ Tqo, 
and only below T/vi there is a weak reduction of (r^), 
in agreement with experiment 22]. In both cases the 
G-OO parameter {t^)g is reduced by singlet orbital fluc- 
tuations, being (t^)g ^ 0.32 (0.37) for LaVOg (SmVOg). 

As a result of the competition between the JT term and 
the CF and orbital-lattice interaction, the temperature 
Too increases first only moderately with decreasing 
and next decreases, resulting in two distinct regimes of 
the phase diagram of Fig. [TJ First, when 8 decreases 
from 157.4° in LaVOg to 144.8° in YVO3, increasing Vat 
dominates and Tqo increases (Fig. [1]). This is similar 
to the i?Mn03 manganites ^ and can be understood 
by considering the total orbital interactions KabTfr^^ in 
the ab planes, including both the superexchange J and 
the JT term Vab, see Fig. HI Here we use again the 
ionic radius vn as in Fig. [1] — we deduced its relation 
to the tilting angle z9, r/j = tq — asin^ 2'd) with tq — 1.5 
A and a = 0.95 A, from the structural data of Refs. 

30, 



21I I22I I . Note that the CF splitting E^ increases 



with decreasing r^j, so it partly compensates the effect 
of increasing Vab- Second, the rapidly increasing orbital 
polarization (Fig. [2]) suppresses the tendency towards 
G-00 and triggers the observed drop of Too (Fig- HI) 
when rn decreases beyond ~ 1.18 A found in YVO3. 
The changes in orbital correlations caused by the 



FIG. 3: (color online) Spin {S^} (solid) and G-type or- 
bital {t^}g (dashed) order parameters, vanishing at Too and 
Tni, respectively, and the transverse orbital polarization (r^) 
(dashed-dotted lines) for LaVOs and SmVOs (thin and heavy 
lines) for Vc = 0.26J; other parameters as in Fig. [2] 



lattice-induced increase of the total orbital interactions 
Kab with decreasing rji (see Fig. H]) suppress the mag- 
netic interactions in the G-AF phase, so the total magnon 
energy scale Wq-af = 4.{Jab + \Jc\) (at T = 0) [9] is re- 
duced from - 1.84J in LaVOa to - 1.05J in YVO3, 
i.e., by a factor close to 1.8. This explains qualitatively 
the rather small magnon energies observed in the G-AF 
phase of YVO3 5]. The reduction is at first instance sur- 
prising as the value of J does not change at all, and it is 
caused solely by the suppression of the singlet orbital cor- 
relations {Ti -fj) by the transverse field Qesi^) (while this 
effect is small for g^ff — 0, in conflict with experiment). 

The role played by the orbital-strain coupling ^ in the 
phase diagram of the i?V03 compounds becomes even 
more transparent by comparing the dependence of g^g 
on the ionic radius with the actual lattice distortion 
u in RVO3 (Fig. [5]). Surprisingly, we find that the ex- 
perimental data for the distortion at zero temperature 
(mo) and above the first magnetic transition (ui) exhibit 
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FIG. 4: (color online) The width of magnon band Wc-af for 
finite Qcs (circles) and without orbital-strain coupling (gcff = 
0, dashed), and orbital interactions in ab planes Kab (squares) 
in the C-AF phase of cubic vanadates (the points correspond 
to the KVO3 compounds of Fig. [T|. Parameters as in Fig. [S] 
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recent observations of coexistence of the G-AF and C-AF 
order [21, 22| will have to be addressed in future work. 



Summarizing, we have introduced a microscopic spin- 
orbital model that provides a satisfactory description of 
the systematic trends for both orbital and magnetic tran- 
sitions in the i?V03 perovskitcs, including the nonmono- 
tonic behavior of the 00 temperature Too- Thereby the 
orthorhombic lattice distortion u, which increases from 
La to Y by one order of magnitude, plays a crucial role 
— it modifies orbital fluctuations and in this way tunes 
the onset of both orbital and spin order in the cubic vana- 
dates. 



FIG. 5: (color online) Experimental distortion (in percent) at 
To = {uo, circles) and above Tni (ui, triangles) for LaVOs 
[20I I and other RYO3 compounds [Tg . [23 | , compared with the 
orbital polarization {r^)T=o and with pcff (O; gcH and g^/K 
are in units of J. Squares show the upper bound for g^/K 
predicted by the theory (at uq — 0). Parameters as in Fig. (2] 
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nearly the same dependence on as either the orbital 
polarization (r^), or the effective interaction ges- This 
is an unexpected outcome of the present theory as infor- 
mation about the actual lattice distortions has not been 
used in constructing the spin-orbital model ([T]), and im- 
plies that: (i) the full i?-dependence of ges is due to u(t?) 
and the bare coupling parameters {g, K} are nearly con- 
stant and independent of rji (Fig. [5]); (ii) g ~ 33 J ob- 
tained using u ~ OmOg^s/J (i.e., g/ao ~ 0.15 eV/A for 
Oq = 3.8 A); (iii) (r^) = x(i^; T)5off(i^), where the suscep- 
tibility X — 0-'2/ J hardly depends on i5 and only weakly 
on T (cf. Fig. 121), so that m(?9) ~ Mo(t?)/[l - xiT)g^/K], 
which justifies a posteriori our approach with a single 
parameter gcs depending only on d; (iv) K > 220 J 
(as X9^/K < !)• K niciy be estimated from the shear 
modulus which is however unknown for i^VOs. Taking 
the data for SrTiOs ^ instead would imply K ~ 6600 J 
(~ 8 eV/A^), i.e., a 3 — 5 % contribution of (t^) to u in 

Soff- 

Finally, we emphasize that the experimental data of 
Fig. [1] are reproduced with rather realistic parameters 
— taking J = 202 K one finds Tmi = 0.73 J = 147 K for 
LaVOg {T"^^ = 143 K 0). Although the present the- 
ory brings us closer to the ultimate understanding of the 
complex phase diagram of the vanadates, several issues 
remain open. One of them is the second phase transition 
at Tn2 to the G-AF phase, which occurs for small rR 0- 
As shown in Ref . Q , this transition is due to an interplay 
between superexchange orbital fluctuations and orbital- 
lattice interactions. While this physics is contained in 
the model ([T]), its quantitative description including the 
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